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ABSTRACT 

LS 5039 is one of a handful of X-ray binaries that have been recently detected at high-energy 7-rays, in this 
case, by the High-Energy Stereoscopy Array (H.E.S.S.)- The nature of this system is unknown: both a black 
hole and a pulsar have been invoked as possible compact object companions. Here we work with a model of 
the high energy phenomenology of the system in which it is assumed that the companion object is a pulsar 
rotating around an 06.5V star in the ^ 3.9 days orbit. The model assumes two different sets of power-law 
spectral parameters of the interacting primary leptons corresponding to the two orbital phase intervals defined 
by H.E.S.S. as having different gamma-ray spectra and very-high-energy (VHE) cutoffs. We show the H.E.S.S. 
phenomenology is completely explained by this model. We present predictions for photons with lower energies 
(for E > I GeV), subject to test in the forthcoming months with the GLAST satellite. We find that GLAST is 
able to judge on this model within one year. 

Subject headings: X-ray binaries (individual LS 5039), 7-rays: observations, 7-rays: theory 



6 



> 

00 



o 
00 
O 



1. INTRODUCTION 

The Gamma-ray Large Area Space Telescope (GLAST) is 
the next generation 7-ray observatory due for launch in May 
2008. Its primary instrument is the Large Area Telescope 
(LAT), which will measure 7-ray flux and spectra from 100 
MeV to -100 GeV. LAT is the successor of the EGRET ex- 
periment that flew onboard the Compton satellite. LAT will 
have better angular resolution, greater effective area, wider 
field of view, and broader energy coverage than EGRET. De- 
veloping predictions for the GLAST energy domain upon 
which to test models of LS 5039 is then timely and worth 
exploring. 

In the previous work (Sierpowska-Bartosik & Torres 2007, 
where we refer the reader for references), under the assump- 
tion that LS 5039 is composed by a pulsar rotating around an 
06.5V star in the ^ 3.9 days orbit, we presented the results 
of a theoretical modeling of the high energy phenomenology 
observed by H.E.S.S.^ The main difference between the pre- 
vious model and the current one is that the former assumed 
a constant spectrum describing the interacting particles along 
the orbit, whereas now it is let to vary. The previous model 
(including detailed account of the system geometry, Klein- 
Nishina inverse Compton (IC), 77 absorption, and cascading) 
was able to describe reasonably well the rich details found in 
the system, both flux and spectrum-wise. However, almost the 
complete spectral energy distribution at superior conjunction, 
as well as other features of the observed orbital variation of 
the H.E.S.S. energy spectra, remained to be consistently ex- 
plained. With models that do not entirely match the observed 
data at the highest energies, the study of lower energy predic- 
tions lacks testing power (strictly speaking, these models are 
already ruled out by H.E.S.S. observations). Here, we ana- 
lyze and motivate the description of the VHE data allowing 
for a variable spectrum of interacting primary leptons along 
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timescale as determined by Casares et al. (2005)- and fast variability dis- 
played on top of this peiiodic behavior, both in flux and spectrum (Aharonian 
et al. 2006). 



the orbit and indicate that the new model can nicely explain 
the H.E.S.S. observations. We then present predictions for 
the flux and spectrum as function of phase for photons with 
lower energies, where H.E.S.S. is not sensitive but GLAST is. 
LS 5039 is near the galactic center and has two other 7-ray 
sources nearby, as well as a significant galactic diffuse back- 
ground: it is not an easy target for GLAST (see Dubois 2006). 
In order to study orbital variability, the need to integrate long 
time intervals and make harder energy cuts to bring the signal 
above background was reported. 

2. THE PULSAR MODEL WITH ORBITAL VARIABILITY IN THE 
INTERACTING PRIMARY SPECTRUM 

The main features of the model are described by 
Sierpowska-Bartosik & Torres (2007). For LS 5039's param- 
eters the star wind dominates over the putative pulsar's, and a 
shock wraps around the latter We assume then that the vol- 
ume of the system is shock-separated as a result of the colli- 
sion between the pulsar and the massive star winds. Between 
the pulsar and the shock, inside the pulsar wind zone (PWZ), 
relativistic leptons are assumed to be frozen in the magnetized 
pulsar wind which propagates radially from the pulsar In this 
PWZ, we compute Klein-Nishina IC cascades against thermal 
radiation from the massive star. Secondary 7-rays move into 
the massive star wind region and whereas some escape the bi- 
nary system, others get absorbed due to 77 process. These 
cascades are followed by means of a Monte Carlo procedure. 
We have computed geometric dependences upon the opaci- 
ties to these processes. For close binaries, the radiation field 
of hot massive stars (type O, Be or WR, having typical surface 
temperatures in the range - IQ-* - 10^ K and linear dimen- 
sion Rs — lOi?0) dominates along the whole orbit over other 
possible fields (e.g., the magnetic field or the thermal field of 
the neutron star). This thermal radiation field is anisotropic 
for e+e^ injected close to the pulsar (the radiation source is 
misplaced with respect to the electron injection place). 

We assume that the injected power in relativistic leptons, 
which initiate the cascading processes, is a fraction of the 
spin-down luminosity, Lsd- We particularize on the sim- 
plest case in which leptons are described, after being repro- 
cessed by losses which dominance can be a function of or- 
bital phase, by a power-law in energy that may be constant (as 
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in Sierpowska-Bartosik & Torres 2007) or vary along the or- 
bit. Here, we focus on the case in which two different power- 
laws are assumed for the interacting lepton population, cor- 
responding to the two orbital intervals proposed by H.E.S.S. 
around inferior (0.45 < 4> < 0.90) and superior ((f> < 0.45 and 
(/) > 0.90) conjunction (INFC and SUPC). Table 1 shows the 
model parameters chosen. In both cases we assume a nomi- 
nal value for Lsd — 10'^^ erg s~^, and i?max=50 TeV. Note 
that LsB and the fraction of it that is converted into relativistic 
leptons (/?) are obviously related. The position of the shock 
(which defines the size of the PWZ where we compute the 
cascades) is defined by the parameter rj — Ls-d/{MVv,c), 
where the wind velocity, V^,, depends on the radial distance 
from the massive star, and M is the star mass-loss rate. Then, 
MVw and Lsd are also connected. Different sets of parame- 
ters (e.g., a smaller Lsd with a higher (3) give similar results. 
This is the consequence of a mild dependence with 77 of the 
distance from the pulsar to the shock. In the models presented 
here, only a small fraction (^1%) of the pulsar's Lsd ends 
up in relativistic leptons. This is consistent with ions carrying 
much of the wind luminosity. 

Once the pulsar injects relativistic leptons (what could itself 
be subject to orbital variability), or alternatively, a close-to- 
the-pulsar shock accelerates a primary population assumed to 
be isotropic in the pulsar rest frame (e.g., Kirket al. 1999), the 
interacting lepton population arises from the equilibrium be- 
tween injection and the losses operating in the system along 
the orbital evolution. In this sense, SUPC and INFC of LS 
5039 may indeed represent qualitatively different physical 
scenarios, as noted already by Aharonian et al. (2006). Let's 
consider first the maximum energy to which leptons are ac- 
celerated in the case the latter proceeds in a shock. The high 
temperature of the star and the hard spectrum of the VHE radi- 
ation measured imply (unless a very hard injection proceeds) 
that IC must happen in the Klein-Nishina regime. For dom- 
inant IC cooling, shock acceleration and cooling timescales 
equality implies that L^max oc (B/w)^'^, where B is the mag- 
netic field and w is the target photon energy density. Since 
w oc d~^, with d the binary separation, and _B is a decreasing 
function of d, e.g., B oc d^^, Ea^^.^ increases by a factor ^ 10 
from periastron (in the SUPC broad-orbital range) to apas- 
tron (in the INFC broad-orbital range), leading to a spectral 
hardening around INFC. At the highest energies, synchrotron 
losses (which timescale is oc L'"^) are expected to dominate 
and steepen the spectrum. The changeover energy from IC to 
synchrotron domination of the radiative losses depend on B 
and d; if B is not exactly oc 1/cZ, it may also introduce a spec- 
tral hardening at apastron. Finally, adiabatic losses happening 
with the typical scales of the system associated with the flow 
patterns -i.e., the stand-off distance of the pulsar wind shock 
(which scales with, and is less than the distance between the 
stars), divided by the post-shock speed (s^ c/3)- generates 
a timescale that is independent of energy and, for LS 5039, 
close to or even smaller than the Klein-Nishina cooling times. 
In general, if leptons enter the production region with a rate 
described by a power-law with index ao, in the case of dom- 
inant radiative losses either by synchrotron or Thompson IC, 
the slope of the interacting population is steppened to + 1- 
In the case of dominance by Klein-Nishina IC, the decrease in 
the cross section is such that the distribution index is slightly 
modified or even hardened at most by < 0.5 in index, with 
hardening decreasing with the ratio Ui^/Us (Moderski et al. 
2005). In the case of adiabatic dominance, the slope of the in- 



teracting population remains the same. Given that the losses 
dominance along LS 5039's orbit can be a function of phase, 
and in addition that also the opacity to pair production is 
phase dependent, it seems natural to assume that in approx- 
imating the distribution of the interacting lepton population 
with a power-law, it has a different index along the two broad- 
phase intervals which qualitative observed features differ the 
most. In this context, the distinction between electron spec- 
tra at INFC and SUPC should be understood as an average of 
a smoother phase-dependence of the electron primary spec- 
trum, what could in turn be tested with future quality of data 
if more complex modeling is available. 

3. MATCHING OF H.E.S.S. RESULTS AND PREDICTIONS FOR 
GLAST 

Fig. 1 shows the results of our model, both flux and 
spectrum-wise, compared with corresponding H.E.S.S. data 
(obtained from Aharonian et al. 2006). The spectra for LS 
5039 is presented in two broad orbital phase intervals around 
INFC and SUPC. Each of the experimental data points repre- 
sents a typical time span of 28 min. The corresponding phases 
for apastron, periastron, INFC, and SUPC are marked. Two 
periods are shown. To ease the comparison with our earlier 
result using a constant spectrum of interacting particles along 
the orbit (Sierpowska-Bartosik & Torres 2007) we include it 
(with a finer gridding) also in these plots. 

H.E.S.S. has also provided the evolution of the normaliza- 
tion and slope of a power-law fit to the 0.2-5 TeV data in 0.1 
phase-binning along the orbit (Aharonian et al. 2006). The 
use of a power-law fit was limited by low statistics in such 
shorter sub-orbital intervals, i.e., higher-order functional fit- 
tings such as a power-law with exponential cutoff were re- 
ported to provide a no better fit and were not justified. To 
directly compare with these results, we have applied the same 
approach to treat the model predictions, i.e., we fit a power- 
law in the same energy range and phase binning. This com- 
parison is shown in Fig. |2] We find a rather good agreement 
between model predictions and data. 

Fig. 1 also shows the spectral energy distribution predic- 
tions extended for energies above 1 GeV. The steeper the pri- 
mary spectrum is, the higher the flux produced at lower ener- 
gies, what is particularly notable for the SUPC broad phase- 
interval (see bold and thin lines of Fig. 1). The minimum flux 
needed for a source to be detected by GLAST after a 1 -month 
and 1-year of operation in all-sky survey, for a E^^ source.** 
At these flux levels, a 20%-uncertainty in the determination of 
the flux, a resulting significance about Su, and a spectral index 
determined to about 6% would be achieved. Even when these 
sensitivities maybe slightly worse for a low-latitude source, 
if this model is correct, GLAST should be able to quantify 
the orbital variability after a few months of integration. The 
fact that the system periodicity is ^^3.9 days allows for a fast 
build-up of integration time around each of the portions of 
the orbit. A month-integration around SUPC will be obtained 
after ~2 months of all-sky survey, putting this model to the 
test soon after GLAST launch.^ Our model also predicts a 
clear anti-correlation between the expected output at TeV and 
GeV energies, which can be seen comparing the lightcurves 
shown in Figs. 1 and 3. The lower the energy range, the 

* www-glast.slac.stanford.edu/software/IS/gIast_lat_perfonnance.htm 
5 LS 5039 is positionally coincident with 3EG J1824-1514 (Hartman et 
al. 1999, Paredes et al. 2000), whose average 7-ray flux above 100 MeV is 
~ 3.5 X 10^^ photons cm^^ s^^. Our model predictions are consistent 
with this observation too. 
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more anti-correlated the signal is, what is a consequence of 
the phase-dependence of the cascading and absorption pro- 
cesses. A hardness ratio defined within the GLAST energy 
domain does not vary as much as it does when constructed 
combining the fluxes at low and high 7-ray energies (a factor 
of 3 versus an order of magnitude), but it may be useful as a 
first check before integration time is achieved for determining 
the spectrum. This is shown in Fig. 4. 

4. CONCLUDING REMARKS 

A detailed modeling of LS 5039, under the assumption that 
it contains a pulsar, is able to fully describe the challenging 
H.E.S.S. -observed phenomenology found in the system. Ob- 
served Ughtcurve, spectra (Fig. 1), and short-timescale spec- 
tral variability (Fig. 2) are matched with a variable spec- 
trum of primary particles along the orbit, where two phase 
intervals are considered. It is interesting to see that the in- 
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fluence of the orbital inclination angle is minor along all en- 
ergy intervals within our theoretical model. Predictions of the 
model at lower 7-ray energy domains are ready for testing 
with GLAST. We find that GLAST is able to judge on this 
model within one year. Features both at the Ughtcurve and 
spectral level can be recognized. Additional tests will come 
at higher energy yet beyond the current reach, e.g., behav- 
ior of the spectra at shorter phase intervals to be compared 
with data from the next generation of ground-based instru- 
ments: H.E.S.S. II, and the planned Cerenkov Telescope Ar- 
ray (CTA). 
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TABLE 1 
Model Parameters 



Parameter 


Adopted value 


Constant lepton spectrum alonfj the orbit 


Fraction of LsD in Icptons: J 
Slope of tlic power-law: i",. 


iO - 
-2.0 


Variable lepton spectrum along the orbit 


Fraction of LsD in leptons at ESTFC interval; (3 
Slope of the power-law at INFC interval: Te 
Fraction of Lgo in leptons at SUPC interval: /3 
Slope of the power-law at SUPC interval: Fe 


8.0 X 10-^ 
-1.9 

2.4 X 10-2 
-2.4 
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Fig. 1 . — Top: Lightcurves for photons with energy above 1 TeV for different binary inclinations (dot-dashed: i = 30 , solid Hnes: i = 60 ), compared with 
H.E.S.S. data. Black lines stand for results obtained with a variable interacting spectrum along the orbit; green (light) lines correspond to the constant spectrum 
case. Bottom: Spectra around INFC (in red) and SUPC (in blue). H.E.S.S. data is shown in the same colors. Results for both cases of interacting electron spectra 
are given. Shaded regions are energy ranges for which we study the lightcurves below. The two horizontal lines between 1 and 100 GeV represent sensitivities of 
GLAST in the all-sky survey mode. 
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Fig. 2. — Shaded ai'eas in the left (right) panel show the change in the normalization (photon index) of a power-law photon spectra fitted to the theoretical 
prediction for each of the 0. 1 bins of phase. The presented results are for the model of the variable interacting spectrum. The two different colors of the shading 
stand for the two inclination angles considered. The size of the shading gives account of the error in the fitting parameters. Data points represent the H.E.S.S. 
results for the equal procedure: a power-law fit to the observational spectra obtained in the same phase binning. 
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Fig. 3. — Predicted theoretical lightcurves for the energy intervals 100 GeV - I TeV, 10 - 100 GeV, and I-IO GeV. Both inchnation angles and interacting 
lepton spectra considered are shown (dot-dashed: i = 30", solid lines: i = 60", black lines: variable spectrum of primary leptons, green (light) lines: constant 
spectrum along the orbit). 



8 



A. Sierpowska-Bartosik & D. F. Torres 




